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ABSTRACT 

As  a  starting  point  for  the  research  about  the 
:'    ~-  on  of  springing,  a  careful  survey  of  the  existing 
literature  on  this  subject  was  made. 

A   mathematical  model  for  the  calculation  of  the  bend.' 

to  springii   i   proposed.   This  mathematical  modi  I 
3  the  modal  analysis  of  the  vibration  of  the  ship  and 
feet  of  the  distortion  of  the  exciting 
.  hu:    or  bhe  calculation  of  the  loading  fi   bion. 

The  influence  of  some  design  parameters  on  the  occurrence 
of  springing  is  pointed  out,  and  it  is  suggested  some  further 

on  .he  pre'"   ,.  with  the  aim  of  evaluating  the 
accuracy  dev  loped. 
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NOMENCLATURE 

A         -  section  area 

A    .  .  .A    -  numerical  constants 
o     n 

A..-,       -  sectional  vertical  added  mass 
B         -  local  beam 

c         -  structural  damping  coefficient  per  unit  length 
(function  of  x) 

-  also  wave  celerity 

C .         - 

1 

C  -   local  sect 

s 

c  -st: 

E  -  modulus  of  elasticity 

EI  -  bending  flexural  rigidity 

G  -sheari 

-  i 

H  -  local  secti    :  aft 

h  -  mean  draft  of  the  sec' 

I  -  sectional  a;:^a  moment  of  ir 

i  -  mode  of  vibration 

I         -  :  length 

r 

I  ( < 

I         -  Bessel  function  of  0  oi 
o 

j         -  imaginary  unit 

K         -  ratio  of  the  a1  '        hear 

at  the  nc  L  L 


k  -  wave  number 

KAG  -  vertical  shear  rigidity 

K.  -  generalized  spring  constant 

L  -  ship's  length 

M  -  bending  moment 

M.  -  weighting  function  for  bending  moment 

M  -  bending  moment  due  to  springing 

N  (x)  -  local  damping  force  coefficient 

N  -  Newman  function  of  0  order 
o 

P(x,t)  -  total  force  per  unit  length  due  to  ship  wave 

interaction 

th 

Q. (t)  -  weighting  forcing  function  for  the  1   mode 


'i 


q. (t)      -  time-varying  beam  deflection  for  the  i   mode, 


in  a  generalized  coordinate  system 
t         --  time 

V  -  ship's  forward  speed 

V  -  shear  force 
s 

V  .        -  weighting  function  for  shear  force 
X. (x)      -  normal  node  shape  for  the  i   mode 

x         -  horizontal  axis  in  direction  of  forward  motion 
of  the  ship  (along  length  of  ship) 

y         -  horizontal  axis  perpendicular  to  x 

z         -  vertical  elastic  deflection 
-  also  vertical  axis 

a         -  angle  between  the  ship ' s  heading  and  the  wave 
direction,  a=0  for  head  seas 


<J>         -  wave  potential 

Y         -  elastic  deformation  angle 

X  -   wave  length 

X  -  effective  wave  length 

X.  -  parameter  defined  by  equation  (28) 

y         -  effective  mass  per  unit  length,  which  is  the 

sum  of  the  ship's  mass  and  the  added  mass  (high 

frequency  limit) 

-  also  coefficient  defined  by  equation  (41) 

y .         -  generalized  mass 

y,         -  frictional  coefficient 

H         -  undistorted  wave  profile 

n*        -  effective  wave  profile 

n         -  unit  vector  normal  to  the  hull  s  surface 

ft         -  Lomrael  function  of  0  order 
o 

(X)  -  circular  frequency 

a)         -  frequency  of  encounter 

th 
to.         -  natural  frequency  of  the  l   mode 

p         -  specific  gravity 

v         -  viscosity 


Dt 

c;i" 
3 


,-)       -  substantial  derivative 
,,-)       -  partial  derivative  with  respect  to  time 


s— )  -  partial  derivative  with  respect  to  space 

)  -  derivative  with  respect  to  time 

) '  -  derivative  with  respect  to  space 

]  -  references  in  Chapter  4 

)  -  equation's  number 

7 


1.   INTRODUCTION 

1.1-WAVE  INDUCED  STRUCTURAL  LOADS 

In  order  to  calculate  the  ship's  hull  structural  res- 
ponse in  a  seaway,  it  is  necessary  to  evaluate  the  loads 
both  in  magnitude  and  in  phase,  so  as  to  make  possible  to 
superimpose  them. 

The  above  mentioned  loads  can  be  generally  considered  in 
two  groups:   the  wave  loads  that  vary  slowly  with  time  in- 
clude the  Stillwater  bending  and  shear,  the  wave  bending 
moment  and  shear,  and  torsion;  and  the  wave-induced  vibratory 
loads,  of  a  higher  frequency,  corresponding  to  the  dynamic 
loading  of  the  hull . 

The  Stillwater  bending  moment  and  shear  is  readily 
calculated  by  available  methods  [1  ] .   The  evaluation  of  the 
wave  bending  moment  and  shear  involves  the  calculation  of  the 
ship's  response  to  regular  and  irregular  seas,  as  well  as 
both  the  short-term  and  long-term  responses. 

When  a  ship  is  moving  in  a  seaway,  five  different  types 
of  forces  contribute  to  the  total  loading:   hydrodynamic , 
hydrostatic,  gravity,  inertia  and  viscous  forces.   It  is 
usual  to  neglect  the  viscous  forces  since  they  are  relatively 
small.   The  hydrostatic,  gravity  and  inertia  forces  havt,  been 
calculated  to  a  good  degree  of  accuracy  for  the  last  few  years. 
The  hydrodynamic  loads  are  of  a  very  complex  nature,  and  their 
evaluation  requires  the  calculation  of  the  pressure  distribu- 
tion around  the  ship's  hull  in  a  seaway.   This  pressure 
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distribution  is  then  used  for  determining  vertical  and 
lateral  bending  moments. 

Theoretical  techniques  [2,3]  are  now  available  to 
evaluate  motions  and  loads  of  conventional  ships  for  five 
degrees  of  freedom  (roll,  pitch,  yavz,  heave  and  sway). 

The  low  frequency  wave-induced  loads  are  considered  in 
the  papers  mentioned  above  in  both  vertical  and  lateral 
directions,  with  the  ship  going  in  any  heading  relative  to 
the  wave  direction.   Lateral  bending  moment,  vertical  bending 
moment,  torsional  moment  and  shear  force  are  calculated  using 
the  results  of  the  motions  of  heave  and  pitch  for  the  vertical 
bending  moment  and  shear,  and  the  motions  of  sway,  yaw  and 
roll  for  the  lateral  bending  moment  and  torsional  moment. 

Torsional  strength  has  become  very  important  for  the 
design  of  containerships  because  of  their  large  hatch  openings, 
therefore  susceptible  to  structural  damage  due  to  torsion 
especially  in  quartering  seaways.   The  analysis  carried  out 
in  [2]  show  the  distribution  of  torsional  moment  over  the 
ship  length.   It  is  important  to  note  that  the  maximum  tor- 
sional moment  doesn't  necessarily  occur  amidships,  but  may 
occur  at  a  considerable  distance  aft  of  it.   There  is  also  an 
important  correlation  between  the  rolling  motion  and  torsional 
moments,  making  the  accurate  calculation  of  the  roll  motion 
a  must. 

In  addition  to  the  low  frequency  wave  loads  discussed 
above,  higher  frequency  wave  induced  vibratory  loads  are  often 
encountered  by  ships  in  a  seaway.   These  vibratory  loads  are 
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caused  by  a  mechanism  totally  different  from  the  one  that 
generates  the  low  frequency  wave  loads,  therefore  resulting 
in  a  completely  independent  response  whose  forces  may,  under 
certain  conditions,  increase  the  maximum  vertical  bending 
moment  to  considerably  high  values. 

This  vibratory  response  is  sometimes  due  to  large  ship 
motions  resulting  in  the  e^^srsion  of  the  bow  or  stern  regions, 
which  leads  to  an  impact  force  when  the  hull  enters  the  water 
(slamming) . 

In  ships  with  large  bow  flare,  the  forces  resulting  from 
the  bow  flare  share  variations  and  the  motions  of  the  ship's 
hull  often  lead  to  a  vibration  known  as  "whipping" . 

We  have  therefore  a  system  of  forces  dependent  on  hydro- 
dynamic  non-linear  effects  and  the  motions  of  a  ship's  hull 
considered  rigid,  with  a  response  that  can  be  approximated 
as  linear.   The  resulting  response  depends  also  on  the 
different  frequencies  and  shapes  of  the  vibratory  behavior  of 
the  hull  (structural  modes) ,  and  appears  as  a  non-stationary 
record  of  high  frequency  oscillations  resulting  from  the 
impulsive  local  forces  at  the  points  of  impact  between  the 
hull  and  the  water  surface.   The  frequency  of  this  vibratory 
motion  is  usually  that  of  the  first  structural  mode  (two- 
node  vibration) ,  because  the  resonance  at  higher  natural 
modes  is  reduced  due  to  structural  and  hydrodynamic  damping. 

It  is  necessary  to  consider  now  the  way  the  wave  bending 
and  slamming  loads  are  superimposed.  Slamming  stress  records 
show  an  initial  effect  of  a  large  compressive  stress  on  the 
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deck  amidships,  in  sagging.   Usually  this  is  not  very  import- 
ant because  it  occurs  when  the  wave  bending  moment  changes 
from  hogging  to  sagging  and  is  therefore  very  small.   The 
"whipping"  vibratory  stress  continues  after  the  slam,  and  a 
critical  situation  may  occur  when  it  is  superimposed  on  the 
bottom  plate  with  the  wave  bending  moment  peak  in  sagging. 
The  critical  condition  for  the  deck  occurs  a  half-cycle 
later,  when  the  maximum  hogging  bending  moment  is  super- 
imposed to  the  remaining  whipping  stress.   It  is  evident  that 
these  conditions  become  more  dangerous  with  large  Stillwater 
bending  moments. 

Besides  the  evaluation  of  the  highest  stress  due  to 
slamming  and  whipping,  it  is  also  important  to  calculate  the 
frequency  of  occurrence  of  slamming  and  the  corresponding 
number  of  stress  reversals  expected ,  so  that  considerations 
on  the  fatigue  performance  of  the  hull  can  be  made. 

Another  vibratory  phenomenon  denoted  as  springing  has  a 
somewhat  different  origin  from  slamming.   It  results  from 
small  motions  of  a  ship  in  relatively  short  waves,  when  the 
frequency  of  encountei  with  the  wave  is  comparable  to  the 
lower  structural  modes  of  vibration  of  the  ship.   Springing 
is  specially  important  in  fast  ships  such  as  destroyers  and 
containerships ,  and  in  large  tankers  and  bulk  carriers.   It 
behaves  like  a  stochastic  process  similar  to  wave  bending 
except  that  it  occurs  at  higher  frequencies.   Springing  can 
be  analysed  as  a  linear  process,  since  it  results  from  a 
resonant  condition  between  the  two-node  vibratory  bending 
response  of  the  ship  structure  and  the  existing  wave. 
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This  phenomenon  will  be  discussed  in  greater  detail  in 
the  next  section. 

1.2-SPRINGING  PHENOMENON 

The  occurrence  of  springing  of  a  ship  is  considered  to  be 
due  to  the  selective  resonance  of  the  ship  hull  vibration  with 
the  frequency  of  encounter  of  the  wave.   It  is  originated  from 
the  synchronism  between  the  hull's  natural  frequency  of 
vibration  and  a  region  of  the  wave  energy  spectrum  of 
sufficiently  high  energy  content  to  excite  a  response. 

As  said  before,  the  particular  ships  for  which  springing 
a  sumes  a  greater  importance  are  large  tankers  and  bulk 
carriers,  as  v/ell  as  fast  ships  such  as  destroyers  and 
containerships .  ■ 

The  evaluation  of  the  wave  induced  vibrations  requires 
the  calculation  of  the  forces  due  to  the  exciting  high- 
frequency  waves,  distributed  along  the  length  of  the  ship, 
and  the  elastic  properties  of  the  ship's  hull  including 
structural  damping. 

The  vertical  hull  force  per  unit  length  is  composed  of 
three  parts:   hydrostatic,  damping  and  inertia  forces.   If 
the  damping  assumes  negative  values,  then  the  vibration  will 
be  self-exciting.   Calculations  carried  out  for  several 
tankers  show  that  the  damping  could  be  negative,  resulting  in 
a  ship  form  more  sensitive  to  springing.   An  extended 
accumulation  of  data  from  the  springing  effect  in  various 
ships  could  eventually  show  a  possible  correlation  between 
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damping  values  and  springing  response.   If  this  negative 
damping  situation  ever  occurs,  most  likely  the  self-exciting 
situation  will  be  only  of  a  transient  nature  and  unstable. 

It  has  been  found  in  previous  treatments  of  the  problem 
of  springing  that  the  damping  (including  the  structural 
damping  and  the  damping  resultant  from  the  distribution  of 
added  mass  along  the  ship  length)  is  relatively  small,  which 
leads  to  a  very  fine  tuning  in  the  response  characteristics 
of  the  phenomenon. 

Analysed  data  [4]  seems  to  show  that  the  springing 
stresses  manifest  themselves  as  a  rather  narrow  band  Gaussian 
phenomenon.   Using  the  superposition  principle  the  ship 
vibrations  in  irregular  waves  can  be  calculated,  if  the 
frequency  characteristics  and  the  energy  spectrum  of  the.  waves 
are  known.   Since  the  ;  pringing  stresses  in  regular  waves 
appear  to  be  linear  with  the  amplitude  of  the  waves,  the 
response  in  irregular  waves  can  therefore  be  evaluated  by 
superposition  of  the  stresses  resulting  from  the  regular  wave 
components  that  constitute  the  irregular  seaway,  as  indicated 
in  [5]. 

The  presenc. -■  of  strong  springing  has  been  observed  even 
in  a  low  sea  sta.be.   When  the  root  mean  square  values  of  the 
springing  bending  moment  are  examined  for  a  particular  sea 
state  it  is  apparent  its  great  sensitivity  to  the  ship's 
forward  speed  and  heading.   The  principal  reason  for  this 
sensitivity  is  that  when  the  ship  changes  heading  or  speed  it 
encounters  a  large  range  of  wavelength  values  that  can  excite 
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the  two  node  natural  vibration  of  the  hull.   This  sensitivity 
is  also  caused  by  the  oscillatory  variation  of  the  amplitude 
of  the  exciting  force,  when  the  wavelength  is  changed. 

Since  the  importance  of  springing  was  recognized,  several 
reports  of  full  scale  high  frequency  stresses  encountered  on 
ships  at  sea  have  been  published,  in  the  hope  of  clarifying 
the  nature  of  the  phenomenon,  as  well  as  provide  experimental 
grounds  to  check  the  various  mathematical  models  proposed . 

Stress  data  from  measurements  carried  out  on  the  Great 
Lakes  ore  carrier  E.L,  Ryerson  show  that  springing  stress 
variations  could  be  as  high  as  64  percent  of  the  wave  bending 
stress  variations  and  magnitudegof  the  springing  stress  on  the 
order  of  14,000  to  15,000  psi  were  found  [6],   A  typical  stress- 
response  spectrum  indicating  combined  effects  of  wave  induced 
and  springing  stresses  is  given  in  Figure  1.   Tests  on  the 
same  ship  also  revealed  that  a  small  change  in  heading  or 
speed  would  change  considerably  the  magnitude  of  the  springing 
stress,  eliminating  it  completely  in  many  cases. 

Usually  springing  stresses  are  observed  much  more  often 
than  long  period  wave  stresses.   Measurements  taken  on  the 
Ryerson  showed  stresses  of  all  magnitudes,  including  stresses 
nearly  twice  as  large  as  the  maximum  low  frequency  wave  stress. 
As  the  springing  phenomenon  has  a  high  frequency  of  occurrence, 
it  is  most  important  to  realize  how  it  can  affect  the  normal 
lifetime  of  a  ship  if  the  springing  stress  is  high  enough  to 
affect  the  fatigue  properties  of  the  structure. 
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The  springing  stress  is  an  important  third  major  stress 
to  add  to  the  normal  still  water  and  low  frequency  wave  bending 
which  have  been  the  loads  considered  in  traditional  design 
methods.   It  was  observed  that  this  phenomenon  was  the  source 
of  the  highest  dynamic  stresses  in  a  normal  shipping  season 
on  Great  Lakes  vessels. 

1.3 -LITERATURE  SURVEY 

Since  the  importance  of  the  high  frequency  wave  excited 
vibratory  loading  of  the  ship's  hull  was  realized,  a  relatively 
large  number  of  works  have  been  published. 

Earlier  papers  on  the  subject  include  those  of  BELGOVA  [7] 
in  1962  and  MAXIMADJI  [8]  in  1967,  which  deal  with  the  vibra- 
tory response  to  regular  waves. 

MATHEWS  [9] ,  in  1967  extended  his  research  to  the  vibra- 
tory response  in  irregular  waves,  with  the  analysis  of  infor- 
mation obtained  on  the  ''Ontario  Power " .   He  suggested  that 
springing  stresses  were  a  function  of  the  geometry  of  the  hull, 
hull  stiffness,  spe~d  and  sea  state,  and  noted  the  significant 
changes  in  springing  response  when  the  ship's  heading  was 
altered.   His  paper  deals  with  the  effect  of  wave  excited 
vibration  on  ship  longitudinal  strength,  in  particular  on  the 
Great  Lakes  bulk  carriers. 

GOODMAN  [10],  in  1970,  developed  a  theory  for  the  predic- 
tion of  springing  based  on  the  assumption  that  the  phenomenon 
is  the  result  of  the  two-noded  hull  vibration  excited  by 
component  waves  of  the  spectrum  that  correspond  in  frequency 
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to  the  natural  frequency  of  the  two-noded  hull  girder  vibra- 
tion.  He  assumes  that  there  is  always  energy  available  in  the 
sea  spectrum  to  excite  the  vibration  at  the  correct  frequency. 
Since  the  period  of  oscillation  increases  with  the  length  of 
the  ship  and  the  energy  level  available  in  the  sea  spectrum 
for  vibration  will  in  general  increase  for  a  fully  developed 
sea,  the  presence  of  significant  springing  in  very  large  ships 
in  fully  developed  seas  should  be  expected. 

VAN  GUNSTEREN  [11]  in  1970,  derived  a  theory  of  wave 
induced- main  hull  vibrations  in  regular  and  irregular  waves. 
The  vibrating  hull  is  represented  by  a  beam-element  model, 
the  wave  exciting  forces  being  determined  with  a  strip  theory 
for  oblique  waves.   The  beam  model  utilized  has  a  discrete 
mass  and  stiffness  distribution.   The  approach  is  somewhat 
similar  to  Goodman's,  but  here  the  bending  moment  is  explicitly 
represented  in  the  modal  form,  and  also  the  responro  at 
frequencies  lower  than  the  two-node  natural  frequency  are 
introduced. 

KUMAI  and  TASAI  [12],  in  1970,  calculate  the  hydrodynamic 
force  induced  by  waves  and  which  results  in  a  high  frequency 
hull  vibration.   The  calculation  is  based  on  the  measured 
hydrodynamic  pressure  at  the  bow  of  a  76,000  d.w.t.  tanker, 
as  a  function  of  time. 

ROBERTSON,  CLEARLY  and  YAGLE  [6]  published,  in  1971,  a 
report  with  the  results  of  strength  studies  carried  over  a 
period  of  four  years  on  the  Great  Lakes  bulk  ore  carrier 
"Edward  L.  Ryerson".   The  purpose  of  this  investigation  was 
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to  obtain  data  of  full  scale  bending  response  and  the  exciting 
waves,  in  order  to  allow  the  study  of  its  correlation.   Along 
with  other  important  conclusions,  the  authors  found  that  the 
vertical  two-node  longitudinal  vibration  (springing)  was  the 
most  frequent  dynamic  loading  of  the  hull  girder.   Dynamic 
stressing  due  to  springing  was  found  in  30  to  45  percent  of  the 
observations,  and  its  magnitude  was  sufficiently  high  to 
deserve  being  included  in  the  loading  considerations  for  design. 

MILES  [4],  presented  in  1971  a  study  on  the  short-term 
statistical  distribution  of  peaks  which  occur  when  low  frequency 
wave  induced  stresses  and  higher  frequency  springing  stresses 
are  present  simultaneously  in  a  ship.   The  observations  took 
place  in  periods  from  approximately  twenty  minutes  to  at  most 
a  few  hours,  so  that  the  sea  stc.te  could  be  considered  a 
stationary  random  p   -■  ss.   Springing  was  found  to  b?.  a  narrow 
d  Gaussian  phenomenon,  with  stress  peaks  highly  correlated. 
Rice  analysis  of  random  noise  [13]  was  successfully  applied  to 
the  springing  phenomenon. 

LEWIS  and  WHEATON  [14]  (1971)  utilized  the  work  of  Miles 
in  the  study  of  the  combined   wave  and  springing  stresses  of 
the  "Edward  L.  Ryerson",  looking  for  a  rational  manner  of 
combining  the  still  water  bending,  the  low  frequency  wave 
induced  bending  and  springing  in  order  to  arrive  at  reasonable 
strength  requirements.   It  was  found'  that  pure  wave  bending 
stresses  are  always  increased  by  springing,  and  that  the  range 
of  maximum  springing  stresses  is  not  affected  greatly  by  the 
range  of  wave  bending  stress.   The  highest  wave  bending  stress 
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and  the  highest  springing  stress  in  a  twenty  to  thirty  minute 
record  seldom  coincide  in  time  and  therefore  the  highest 
combined  stress  in  any  record  is  alw    less  than  the  sum  of 
the  highest  wave  bending  stress  and     highest  springing 
stress  in  the  record.   The  authors  o:   jferenoe  [14]  point  out 
the  feasibility  of  predicting  the  ma:  Lmum  combined  stress  when 
the  separate  wave  and  springing  stresses  are  known  or  calculated, 
and  from  there  a  long-term  prediction  of  the  combined  stress, 
as  previously  done  for  wave  bending  only. 

KAPLAN  and  SARGENT  [15]  (1972)  evaluated  the  springing 
response  assuming  that  the  bending  moment  (due  to  springing) 
is  represented  in  terms  of  the  hull  acceleration  associated 
with  the  mode  shape,  hence  eliminating  significant  contribu- 
tions at  lower  frequencies.   The  results  allow  the  presentation 
of  a  time  history  of  the  vertical  bending  moment,  for  a  parti- 
cular exciting  wave  system  time  history,  as  well  as  a  repre- 
sentation in  spectral  forms  via  frequency  response  methods. 
The  root  mean  square  values  obtained  using  the  two  methods 
are  in  agreement  for  the  same  wave  system,  and  therefore  the 
choice  of  any  one  of  the  approaches  depends  upon  the  particular 
use  of  the  results  (for  example,  the  spectral  form  will  be 
used  for  long  term  prediction,  and  the  time  history  for 
comparison  with  details  of  experiments  involving  precise 
phase  relations) .   The  solution  of  the  hydrodynamic  problem 
of  high  frequency  local  wave  loading  is  attempted  using  a  two- 
dimensional  section  approach,  but  in  spite  of  its  different 
mathematical  form,  the  numerical  results  obtained  show  very 
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little  difference  from  the  treatment  of  present  strip  theory- 
procedures  for  the  case  of  a  semi-circular  cross  section. 

KUMAI  published  in  1972  a  paper  [16]  on  the  ship's 
hull  vibratory  response  to  high  frequency  wave-induced  loading, 
followed  by  another  work  [17]  in  1973  in  which  the  effect  of 
distortion  of  the  wave  on  the  ship  side  is  included.   In  the 
later  work  the  wave  induced  vibratory  force  and  response  are 
evaluated  using  the  ship, side  wave  height  as  measured  by  self- 
propelled  model  test  of  a  tanker  in  regular  waves.   The  results 
seem  to  show  that  the  exciting  force  is  distributed  along  the 
ship  length,  although  concentrated  on  the  region  of  entrance 
of  the  hull  because  of  the  effect  of  the  distorted  wave,  which 
in  this  case  was  two  to  three  times  higher  at  the  bow  and  a 
half  to  a  quarter  lower  at  the  stern  than  the  undistorted  wave. 

1  :         :  I  hor  also  derives  an     re:   ;.on  that  shows  that  a  ship 
with  large  beam  to  depth  ratio  will  be  subjected  to  v. 
severe  springing  stresses. 

KAPLAN  et  al  [18]  and  NORDENSTROM  et  al  [19]  published  in 
197  3  the  reports  of  the  Ship  Structure  Committee  (Committees 

2  and  3  respectively) ,  where  recent  developments  on  the  evalua- 
tion of  wave  loads  are  dis<       to  a  great  extent.   The  Report 
of  Committee  2  concentrates  on  the  hydrodynamic  aspects  of  the 
wave  loading  (both  low  and  high  frequency) ,  while  the  Report 

of  Committee  3  offers  a  statistical  approach  to  the  problem. 
Both  reports  offer  extended  lists  of  references. 


20 


Tests  on  springing  are  being  carried  out  at  the  Webb 
Institute  of  Naval  Architecture  in  order  to  continue  the 
investigation  using  models  and  allow  further  computations  in 
an  effort  to  obtain  some  correlation  between  theory  and 
experiment.   The  Shipbuilding  Laboratory  at  Delft  is  planning 
tests  with  a  segmented  model  of  a  tanker  in  order  to  obtain 
data  on  springing.   Hopefully  the  results  from  these  and  other 
tests  will  help  the  understanding  of  the  fundamental  parameters 
that  affect  the  high  frequency  vibratory  response  of  ships  and 
also  the  sensitivity  of  the  results  to  the  assumptions  used  in 
current  mathematical  models. 
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2.   MATHEMATICAL  MODEL 

2 . 1-GENERAL 

We  will  proceed  now  to  find  a  system  of  equations  that 
will  provide  an  adequate  description  of  the  ship  structural 
response  to  springing.   The  two  following  steps  would  be  the 
conversion  of  these  equations  into  a  computer  program  and  the 
evaluation  of  the  proposed  mathematical  model  by  comparing 
its  predictions  with  experimental  results;  these  two  steps 
arebeyond  the  scope  of  the  present  work  and  therefore  will 
not  be  carried  out  here. 

When  considering  the  problem  of  the  vibration  of  a  ship 
hull,  it  is  usual  to  consider  the  hull  structure  as  an  elastic 
beam  with  free  ends  and  nonuniform  mass  along  the  ship's 
length.   The  vibration  is  therefore  governed  by  the  Timoshenko 
equation  incorporat       iding,  shear  and  rotary  inertia. 

The  numerical  methods  that  have  been  used  in  the  solution 
of  this  beam  problem,  when  compared  with  data  taken  from 
experiments  aboard  ships,  show  in  general  a  good  agreement, 
except  for  higher  modes  of  vibration  (above  the  fourth)  where 
the  measured  natural  f  requenci'.:  ?.  tend  to  be  lower  than  the 
predicted  ones.   The  .  >  '  tc ;  a.1  reasons  that  have  been  pointi 
out  to  explain  such  discrepancies  are  the  deficient  evaluation 
of  parameters  such  as  stiffness,  load  distribution,  added  mass, 
which  lead  to  inaccurate  calculations,  and  the  possibility 
that,  after  a  certain  mode,  the  ship  does  not  behave  as  a 
beam  anymore  [20] .   Other  possible  reasons  are  the  shear  lag 
effect,  not  taken  into  consideration  by  the  beam  theory,  and 
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the  fact  of  neglecting  the  coupling  between  horizontal,  ver- 
tical and  torsional  motions,  as  well  as  local  vibratory 
responses  [21] . 

Other  beam  models  have  been  experimented  to  study  the 
vibratory  response  of  ship  hull,  such  as  the  three-dimensional 
beam-shell-sprung  body  in  which  the  shafting  and  propeller  is 
modeled  as  a  sprung  body  connected  to  a  beam-shell  representa- 
tion of  the  hull.   Other  examples  are  those  of  Kline  and 
Clough  [22]  who  used  two  beams  interconnected  by  rigid  links 
and  springs  and  the  one  of  Ghtaka  [23] ,  who  used  two  elasti- 
cally  connected  parallel  beams. 

The  finite  element  approach  proposed  by  some  authors  has 
the  advantage  of  representing  the  hull  by  a  three-dimensional 
structure,  therefore  providing  a  model  closer  to  reality, 
without  the  inconveniences  of  the  one  dimensional  method, 
and  automatically  takes  into  account  the  effects  of  shear  lag, 
local  vibratory  motions  and  coupling  with  torsional  vibrations, 
Unfortunately  the  finite  element  method  requires  a  much  more 
elaborate  computation  technique,  because  of  the  need  of  three- 
dimensional  stiffness  factors,  wave  loads  and  added  mass; 
besides,  the  irapsrovement  in  accuracy   for  lowor  modes  of 
vibration   is  small,  and  does  not  seem  to  justify  the  incre;  : 
in  computing  time. 

From  what  was  said  above  it  is  concluded  that  the  beam 
theory  is  a  satisfactory  representation  of  the  vibratory 
response  of  the  ship's  structure,  provided  that  no  significant 
response  is  expected  from  the  higher  modes  of  vibration. 
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2. 2 -MODAL  ANALYSIS 

The  modal  analysis  for  the  study  of  the  bending  moment 
has  the  important  advantage  of  leading  to  a  mathematical 
formulation  not  too  cumbersome,  if  certain  simplifications 
are  assumed.   In  the  modal  analysis  the  equation  variables 
are  expressed  as  a  product  of  two  functions,  one  depending  on 
time  only  and  the  other  on  space  only,  and  the  solution  is 
represented  in  terms  of  a  series  of  normal  modes.   It  is 
required  an  independent  solution  for  the  normal  modes,  which 
are  determined  by  solving  the  eigenvalue  problems  for  the 
natural  frequencies  and  the  mode  shapes  (eigenf unctions) . 

The  vibratory  response  of  a  ship  structure  considered  to 
be  an  elastic  beam  with  non-uniform  mass  and  free  ends  is 
governed  by  the  following  partial  differential  equations, 
which  have  been  already  presented  in  [24,  25] : 

y     +  c    +   a   =  P(x,t)  (1) 

dx 

where : 

t    -  time 

Vo    -  Shear  force  (function  of  x,t) 

x    -  horizontal  axis  in  direction  of  forw.  rd  motion 

of  the  ship  (along  length  of  ship) 
z     -  vertical  elastic  deflection 
c     -  structural  damping  coefficient  per  unit  length 

(function  of  x) 
P  (x,  t) -total,  force  per  unit  length  due  to  ship-wave 
interaction 
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\i  -  effective  mass  per  unit  length,  which  is  the  sum 

of  the  ship's  mass  (per  unit  length)  and  the  added 
mass  (per  unit  length).   Function  of  x  only,  based 
on  the  ship's  still  water  line  (high  frequency 
limit) . 

TT  =   V   +Ir&  (2) 

9x     s     i  3t 

where : 

I    -  mass  moment  of  inertia  per  unit  length  with 

respect  to  an  axis  normal  to  the  xy  plane  (function 

of  x) 
M    -  bending  moment  (function  of  x,t) 
Y    -  elastic  deformation  angle 
Ir(32y/3t2)  -  rotary  inertia 

M  =  EI^  (3) 

where : 

E  -  modulus  of  elasticity 

I  -  sectional  area  moment  of  inertia 

EI  -  bending  flexural  rigidity 

a       V 

Tfx  "  "  KAG   +  Y  (4) 

where: 

KAG   -  vertical  shear  rigidity 

K    -  ratio  of  the  average  shear  stress  to  the  sheer 
stress  at  the  neutral  axis  under  vertical  load 
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A    -  section  area 
G    -  shear  modulus 

This  system  of  partial  differential  equations  has  to 
satisfy  the  boundary  conditions  at  the  free  ends: 

Vg  (-L/2,  t)  -  Vs(L/2,  t)  =  0 
and  (5) 

M(-L/2,  t)  =  M(L/2,  t)  =  0 

Reference  [24]  discusses  different  possible  procedures 
for  the  solution  of  this  system  of  equations  for  the  problem 
of  bow-flare  slamming.   It  was  found  -that  the  methods  of 
direct  solution  consisting  in  converting  the  partial  differ- 
ential equations  into  ordinary  differential-difference  equa- 
tions, by  breaking  the  berim  into  a  large  number  of  modal 
egments,  were  inadequate  because  of  computational  limitat:  ns, 

This  method  of  handling  the  parameters  and  variables  of 
the  system  of  equations  as  discrete  rather  than  distributed 
has  also  the  disadvantage  of  requiring  the  evaluation  of  the 
bending  and  shear  rigidity  distributions,  for  which  an  exact 
procedure  is  not  readily  available.   This,  along  with  the 
approximate  numerical  method.1  of  solution  of  the  system  of 
partial  differential  equations,  becomes  the  source  of 
important  inaccuracies. 

An  alternate  procedure  is  possible  if  the  term  that 
includes  the  rotary  inertia  I   can  be  neglected.   In  this  case 
the  system  of  equations  can  be  reduced  to  a  single  second 
order  linear  differential  equation,  with  constant  coefficients 
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and  where  the  parameters  EI  and  KAG  (which  are  both  functions 
of  x)  are  not  included. 

Since  the  rotary  inertia  has  a  reduced  influence  on  the 
predominant  mode  of  ship  elastic  response,  which  is  the  first 
mode  of  motion,  the  beam  can  be  represented  by  a  modal  model. 
Reference  [20]  reports  experiments  carried  out  on  destroyers, 
showing  that  the  effect  of  the  rotary  inertia  on  the  natural 
frequency  of  the  vertical  vibration  is  reduced.   In  [11]  the 
author  states  that  for  prismatic  bars  the  influence  of  the 
rotary  inertia  on  the  natural  frequency  is  a  quarter  of  the 
influence  of  shear,  this  being  very  small  when  compared  with 
bending,  and  therefore  carries  out  the  analysis  without 
including  the  rotary  inertia. 

Under  the  assumption  that  the  rotary  inertia  can  be 
neglected,  it  is  possible  to  prove  that  the  dynamic  behavior 
of  the  beam  can  be  treated  in  terms  of  a  series  of  re   )nses 
in  each  of  its  normal  modes  i.   Equation  (6)  represents  the 
orthogonality  property  between  the  effective  mass  per  unit 
length  and  the  normal  modes: 

L/2 

/_L/2  \i(x)     '    X±(x)  '  X.  (x)  -  0  (6) 

th 
where  X. (x)  is  the  normal  mode  shape  of  the  i"  mode;  it 

represents,  in  arbitrary  dimensionless  units,  the  shape  of 

the  relative  displacement  along  the  length  of  the  beam  for 

the  mode  i . 

The  different  variables  in  the  system  of  differential 
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equations  are  represented  in  product  form  as 


z(x,t)  =   Z    q. (t)  •  X. (x)  (7! 

i=l    x  x 


M(x,t)  =  Z         q. (t)  ■  M. (x)  (8) 

i=l  x 


V(x,t)  =   Z   q  (t)  ■  V    (x)  (9) 

S         i=l   X        Si 


In  the  preceding  equations  q. (t)  is  the  time-varying 
am  deflection  for  the  i   mode,  in  a  generalized  coordin; 
system. 

Hence,  equation  (7)  represents  the  summation  of  the 
contributions  from  all  the  modes,  obtained  for  a  particular 

mode  i  by  multiplying  the  time-varying  beam  deflection  for 

.  th 

'-  i  1J  mode  q.  (t)  by  the  dime   ionless  nor:  '■ 

X. (x).   Similarly  equation  (S)  represents  the  bending  moment 
M(x,t)  as  the  product  of  q. (t)  by  a  weighting  function  M. (x) 
and  equation  (9)  gives  the  shear  force  V  (x,t)  as  the  product 
of  q. (t)  by  the  weighting  function  V  . (x) .   The  spacial 
functions  M. (x)  and  V  . (x)  will  be  determined  later. 

1  S  .1. 

The  leading  function  can  be  represented  by  the  f<  1] 

serie.'  expression: 

U(x) -Q. (t) -X. (x) 
P(x,t)  =   _Z   -I-72 (10) 

1=1  /_L/2y (x) -X|(x) -dx 


th 
where  Q. (t)  is  the  weighted  forcing  function  for  the  i   mode, 
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and  can  be  found  if  we  multiply  both  sides  of  (10)  by  X. (x) 

and,  using  the  orthogonality  property  (6),  integrate  over  the 

ship' s  length: 

L/2 
Qi(t)  =    f-L/2      p(x't)  '  xi(x)  dx  (11) 

Combining  the  fundamental  elastic  equation  (3)  with  the 
relation  between  bending  and  shear  effects  (4),  neglecting 
the  rotary  inertia  I_,  and  substituting  equations  (7)  through 
(10)  into  the  initial  equations  (1)  through  (4),  the  following 
system    differential  equations  is  obtained: 


VQ.X. 
Z   [n.q.X,  ■!-  cq.X.  +  q,  •  (V,  )  '  -   --^ ]  =  0    (12) 


i=l 


li  i    ^i"i    ^i  ^i'       L/2 


/   /0  yX?dx 
—  L//C.  X 


X: 


I  ■>'.  ■  c.  {:      )   '  j  -  0  (13) 

'  X  X      "X    X 


and 

00  V.       q  .M. 

£   [q.-(X.),!    +  q.  •  (=|tt)  ■ i~]  =  0  (14) 

.  ,   ax    x        Jx     KAG       Ex 

ere  the  dot  means  differentiation  with  respect  to  time  and 
the  prime  means  dif:        bion  with  respect  to  space  (x) . 

Equations  (12),  (13)  and  (14)  can  be  combined  into  the 
following  expression: 

yq.X.  +  cqiXi  +  q.   •   [BI.^)-  +  EI  ■   <J^)  •  ]  »   - 

•  vQ.x. 
s— i-i  (15) 
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where : 

L/2 

m-±  -  /_L/2  yx*  dx  (16) 


Considering  the  free  motion  of  the  beam,  with  no  forcing 
function  present,  equation  (15)  can  be  expressed  as 

qi+(c/iJ)qj_       1  V.        (M.  )" 


±-    [EI-(X,  )"  +  EI-  (_i_)  •]'•=— -±—  (17; 


q.  yX.  L    v  l'         KAG'  J     yX. 


It  will  be  discussed  later  that  c/y  can  be  assumed 
constant;  therefore,  equation  (17)  contains  a  function  of  t 
alone  on  the  left  side,  and  the  right  side  is  a  function  of  :■• 
alone.   Using  the  normal  procedure  of  separation  of  variables 
on  partial  differential  equations  each  side  is  set  equal  to  a 
constant.   In  this  case,  to  assure  stable  responses,  the 
nt  is  -to?,  where  co.  is  the  natural  frequency  of  the 

11  i  J. 

mode . 

Setting   the    first   term   in    (17)    equal    to    -co?    gives 

a.    +    (c/y)q.    +    co?    q.    =    0  (18) 

hich  is  a  linear  second  order  differential  equation  with 
constant  coefficients. 

Setting  the  last  term  in  (17)  equal  to  -co?  loads  to 

(Mi)"  =  y  w*  X±  (19) 

which  after  integration  over  the  ship's  length  gives  an 
expression  for  the  spatial  weighting  function  M. (x) : 
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M.  -  fX  fX  y  03?  X.  dx  dx  (20) 

1     -L/2    -L/2      1   X 


This  double  integral  can  also  be  expressed  as  the  following 
single  integral,  which  is  a  more  convenient  form  for  M. : 

M.  =  w?  x  (x-s)  ■  y(s)«  X.  (s)  ds  (21) 

i     l  j_L/2  i 

Returning  to  the  forced  motion,  we  can  combine  equations 
(17)  and  (19)  into  (15)  to  obtain: 

y.  q.  ■!-  C  q.  +  K.  q.  =0.  (t)  (22) 

where  the  generalized  mass  y.  is  defined  by  equation  (16), 
the  generalised  damping  C,    is 

C.  -  fL/2        c-X?-<  (23) 

-L/2 

and  the  generalized  spring  constant  K.  is 


K.  =  u?*y,  (24) 

ill  v   ' 


From  equation  (16)  we  Fee   that  y,  is  constant  for  a 
particular  mode  i,  because  y,  the  effective  n  tss  per  i  ait 
length,  is  only  a  function  of  x. 

To  define  the  damping  coefficient  associated  with  hull 
vibration  we  will  use  here  the  approach  followed  in  [26] f 
which  considers  three  factors  contributing  to  damping:   a 
structural  damping  factor  per  unit  length  c  /y ,  function  of 
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the  mode  frequency  only,  where  minor  hydro dynamic  effects  of 
water  friction  and  generation  of  pressure  waves  are  also 
included;  the  effect  of  generation  of  surface  waves;  and  a 
correction  factor  due  to  forward  speed,  -V(dA,_/dx),  where 
A_?  is  the  sectional  vertical  added  mass. 

Since  the  ship's  structural  natural  frequencies  are 
relatively  high,  the  effect  of  damping  due  to  wave  generation 
can  be  neglected  when  compared  with  the  other  factors 
included  in  c^/ii.      We  have  therefore  the  following  damping, 
per  unit  length? 

dA 

c  =  c   -  V  —3^-  (25) 

s      dx 

which  when  multiplied  by  the  square  of  the  normal  mode  shape 
of  the  i   mode  X?,  and  integrated  over  the  ship's  length, 


.  ■    to 


s.-    ,r   L/2      3  3    v2 


C.  =  (~)Vi    -   V  ^^  -£L    ■  Xf(x)  dx  (26) 


The  value  that  should  be  used  for  the  generalized  spring 
constant  K.  is  the  one  obtained  if  the  high  frequency 

limit  of  the  added  inuss  is  used  in  (24)  because  it 
is  assumed  that  thi s  value  has  already  been  reached  at  the 
relatively  high  ship's  structural  mode  frequencies. 

From  the  above  discussion,  and  since  the  natural  frequency 
of  the  i  '  mode  co  •  has  a  fixed  value,  we  conclude  that  y  .  ,  C. 
and  K.  are  constants.   Equation  (22)  is  therefore  a  linear 
second  order  differential  equation  with  constant  coefficients 
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and  a  forcing  function  Q. (t) . 

For  the  sake  of  completeness  the  closed  form  solution  of 
this  equation  is  given  in  (27)  for  the  initial  conditions 
q. (0)  =  q. (0)  =  0,  i.e.,  assuming  that  at  t=0  the  ship  is  at 


rest, 


,   Q.  (t)    _c_  ,    . 
q.  (t)  =  r      yi- [e  2p  u  t']-  sin [X     (t-T) ] dx  (27) 

1        o   Aiyi  X 


where 


X.    =   Voj?  -  i-  (%2  (20) 

yi 

and  t  is  the  dummy  inter:      i  variable. 

For  the  springing  phenomenon,  the  solution  for  the 
vibratory  response  in  the  frequency  domain  of  equation  (22) 
is  readily  obtained  if  bhe  forcing  function  Q. (t)  is  rep: 
sented  in  the  following  form: 


Q±(t)  =  Qq  exp(j  oje  t)  (29) 

where  j  is  the  imaginary  unit  and  the  frequency  of  encounter 

u   is,  for  head  seas, 
e  ' 

w   =  ~    (V+c)  (30) 

e    X 

where  X    is  the  wave  length,  V  is  the  ship  forward  velocity 


and  c  is  the  wave  celarity  (  /gX/2iT)  . 

The  steady  state  solution  of  equation  (22)  is  therefore 
given  in  complex  form  by 
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Q~  exp(j  co   t) 

q   =  -°— - S (31) 

K.-ury  .+j  co  C.  - 
1   e  i     e  l 

which  allows  the  determination  of  amplitude  and  phase  as  a 

function  of  the  frequency  of  encounter  co  . 

^    --  e 

Equation  (21) ,  which  was  derived  assuming  that  the  main 

contribution  to  the  bending  moment  comes  from  the  inertial 

loads  along  the  hull  resulting  from  the  vibratory  deflections 

and  including  the  fluid  inertial  force  associated  with  the 

added  mass,  should  nc>\>   be  substituted  by  the  followJ 

expression,-  for  the  midship  bending  n  anient: 

o 
Ms  =  -w?   /_L/2  x-y(x) •Xi(x)-dx-qi(t)  (32) 

The  above  expression  has  only  a  significant  response  on 

irst  mode  of  natural  vibration  of  the  ship  structure 

(for  i-  1) ,  hence : 

o 
Ms  =  f_L/2   x-y(x)-Xi(x)-dx-qi(t)  (33) 

I   ause  the  influence  of  the  inertial  reactions  associate 
with  the  ship  structural  deflection  accelerations,  at 

quencies  higher  bh    .i   first  raode,  will  contribute  to 
an  increase  in  the  bending  moment  response  for  the  first 
mode.   On  the  higher  range,  the  decay  of  the  wave  forces 
exciting  the  hull  deflection  will  greatly  reduce  the  springing 
bending  moment  response.   The  response  should  then  be  given 
for  the  first  mode  by  equation  (33)  ,  with  a  rapid,  asymptotic 
decay  for  higher  natural  modes  of  vibration. 
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2.3-LOADING  FUNCTION 

In  order  to  evaluate  the  wave  induced  bending  moments  it 
is  necessary  to  find,  along  the  ship  length,  the  distribution 
of  the  loads  due  to  local  wave  forces  as  well  as  those  due 
to  the  rigid  body  motions  of  the  ship. 

It  is  assumed  that  the  influence  of  the  ship  rigid  body 
motions  on  the  springing  phenomenon  is  minor,  so  that  the 
loading  force  responsible-  for  the  high  frequency  vibration 
results  only  from  the  interaction  between  the  waves  and  the 
sbJ  ;;■' :  hull « 

The  force  per  unit  length  i«  due  to  hydrostatic  and 
hydrodynamic  effects.   The  hydrodynamic  loading  includes  the 
inertia  force  due  to  body  dynamic  motions  and  the  dissipative 
damping  force.   The  hydrostatic  term  accounts  for  the  exciting 
due  to  the  oncoming  waves. 

The  local  vertical  wave  force  acting  on  a.  section  of  a 
ship  is  then  given  by  [10] . 

P(x,t)=-CpgBn*  +  [N  (x)-V-g~]n*  +  A-  fi*>e  (34) 

where ; 

-  local  beam 

N  (x)  -  local  damping  force  coefficient 

V  -  ship's  forward  speed 

A^_  -  local  section  vertical  added  mass 

n*  -  effective  wave  profile 

h  -  mean  draft  of  the  section,  and  can  be  approximated 


by 
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h  =  HC  (35) 

s 

where : 

H     -  local  section  draft 

C     -  local  section  area  coefficient, 
s 

The  derivatives  r\*   and  r\*    should  be  calculated  as  total 
derivatives,  since  the  rate  of  change  occurs  not  only  in 
time  but  also  in  space:  • 

n*  =  §n!L  .  (3   v  5— >n*  (36) 


*•-&  -   'It  -v  fe^  (37, 


Attention  is  now  focused  on  the  evaluation  of  t; 
fective  wave  profile  ti*. 

Previous  evaluations  of  springing  response  to  wave 
induced  high  frequency  loading  have  been  made  considering 
the  undistorted  oncoming  wave  as  the  water  profile  responsible 
for  the  excitation.   It  is  known  that  when  a  ship  crosses  a 
there  is  a  very  important  distortion  of  the  wave  prof 
the  water  surface  geom  try  at  the  ship  side  being  far  from 
similar  to  the  undisturbed  wave  profile. 

The  effect  of  this  distortion  is  even  more  important  if 
we  notice  that  the  amplitude  of  the  distorted  wave  is  reduced 
as  it  progresses  towards  the  stern  of  the  ship,  therefore 
creating  at  the  bow  a  region  where  the  load  is  concentrated. 
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Realizing  the  importance  of  this  phenomenon,  Kumai  [17] 
carried  out  a  series  of  experiments  with  a  self  propelled 
model  tanker  in  order  to  determine  the  distortion  of  the  water 
line  and  the  profiles  and  envelopes  of  the  wave  at  ship  side. 
Using  this  data  a  loading  function  was  defined,  from  which 
the  acceleration  response  at  the  bow  of  the  ship  was  evaluated 
and  with  the  bow  acceleration  the  stress  amidship  was 
calculated  [16] . 

An  attempt  is  made  here  to  include  a  mathematical  evalua- 
tion of  the  shape  of  1    distorted  wave  at  ship  side  in  the 

ral  procedure  for  the  calci   :  Lor  c   the  springing  res- 
ponse.  The  outline  follows  the  analysis  used  by  GRIM  [27]  in 
1957  for  the  calculation  of  wave  generated  forces  on  a  ship 
hull. 

When  a  ship  cross  s  a  i  ave  it  is  clear      the  forces 
produced  by  the  hull  on  the  water  surface  will  generate  a 
deformation  which  is  similar  to  the  exciting  wave,  oscilla- 
ting both  in  time  and  in  space  (along  the  length  of  the  ship) . 
If  the  wavelength  is  relatively  large  compared  with  the 
length  of  the  ship,  then  there  is  no  qualitatively  large 
ference  from  the  heaving  motion  in  still  water*   The 
ces  produced  by  the  shiphull  on  the  water  surface  will 
therefore  generate  a  system  of  circular  waves  with  a  length 
which  is  about  the  same  as  the  length  of  the  exciting  wave, 
and  because  the  magnitude  and  the  direction  of  the  forces 
produced  by  the  shiphull  are  not  approximately  constant  over 
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the  length  of  the  ship  (as  it  happens  in  the  case  of  the 

heave  motion) ,  but  change  in  the  same  way  as  the  exciting 

wave,  the  amplitude  of  this  exciting  wave  falls  off  along 

the  ship,  resulting  in  a  strong  deformation  in  amplitude, 

while  the  frequency  of  the  oncoming  wave  is  not  significantly 

altered. 

This  deformation  has  to  be  taken  into  account  in  the 

calculation  of  the  forcing  function  (34),  and  to  evaluate 

the  shape  of  the  distorted  wave  at  the  ship  side  the  following 

wave  potential  in  tw  i  dj  nensions  will  be  used  as  a  starting 

point  (the  time  function  e  ' c  is  left  out) : 

m      exp(-Kz)cos(Ky) 
<Ky,z)  =  A  /   — „  .  ,  . dK  + 


CO 

1  ik 

I      A  Re  -{ i-=- — —  • 

n=l  n     (y+iz)      (2n-l) (y+iz)^n~ 


where  k  is  the  wave  number 

e    2tt 
k  =  —  -  ~  (39 

g     A 

and  u,  is  a  frictional  coefficient  [28] 

,J1  =  —  =  'V  gA  (40) 


with 

y  =  v  p  (41) 
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and 

v     -  viscosity  (L2T  1) 

p     -  specific  gravity  (ML  3) 

c     -  wave  celerity 

X  -   wave   length 

g     -  acceleration  of  gravity 

The  frictional  coefficient  y,  is  always  small  [29]  and 

usually  neglected  [30] . 

The  constants  A   in  (38)  are  determined  in  such  a  way 
n 

that  the  conditions  at  the  boundary  are  satisfied,  or 

■  \   =  0  (42) 

on  the  hull  surface,  where  ri  is  the  unit  vector  normal  to  the 
surface. 

In  three-dime]      ,  the  e    bion  .    the  potential  wil! 
be: 


exp(i  w^t)     °°     L/2 

(43) 


4>(x,y,z)  =   = I        I  A(Fj-ib[(Z~x)ry,z]d^ 

*  n=0    -L/2   n     n 


where : 


V 


°     ,/Tc- 


.".-  e  '   {No[k/Tl   )  ..     : 


(5-x) 2+y2+z2 


+  R  [k/(£-x)  2+y*]  +  -  /       eku  du  + 


/ 


2  ,„2  .,,2 


(£-x)  ^H-y^+u 


+  2i  I   [k/(£-x)2+y2  ]1  (44) 
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for  n>o 


r.2  (n-1)        -.   z  -,  , 

^   =  *     _  -[ ±z -j- ilii? —-y    (45) 

9z        [  (t,-x)  +y  +z  ]       [  (q.-x)  +y  +z  ] 


and  [31] 

I      -  Bessel  function  of  0  order 
o 

N     -  Newman  function  of  0  order 
o 

0,  -   Lommel  function  of  0  order 

o 

The  potential  (43),  for  any  function  A  (£)  (C  is  the 
dummy  coordinate  associated  with  the  x  direction),  satisfies 
the  continuity  cc     Lon  ;  d  bhe  free  surface  condition.   T 
constants  A  (£)  are  calculated  from  (38)  so  that  the  boundary 
conditions  on  the  shiphuil  are  satisfied,  and  they  must 
disappear  beyond  the  hull's  length  (j£.  |>L/2). 

me1  >d  is  fc  srefore  based  upon  the  determine 
a  source  distribution  for  the  two-dimensional  cross-s    Lons, 
the  three-dimensional  problem  being  solved  from  the  combina- 
tion of  the  pressure  distribtuion  obtained  in  two-dimensions ; 
the  results  should  constitute  an  improvement  over  the  usv 
strip  method. 

The  potential  of  the  undii  bu    :  >..  re.   must  be: 

<j>  =  -2  exp[iwt  -i-  ik(x+iz)]  (46) 

W    OJ 

If  we  take  off  the  time  factor: 

<b  =  ^  exp[ik(x+iz)]  (47) 

w   OJ 
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The  amplitude  of  the  undisturbed  wave,  using  as  reference  the 
plane  z=0,  will  therefore  be: 

n  =  helkx  (48) 

For  similarity,  in  the  initial  equation  (43)  for  the 
three-dimensional  potential  of  the  distorted  wave  due  to  the 
ship  hull,  we  will  use  instead  of  the  function  A  (£  )  : 

A  {£)+  elk   ■  A  (5)  (49) 

n  ii 

This  can  be  done  because  th  :  conl  Lnuity  condition  ano'  the 
conditions  on  the  f:  surface  are  not  violated.  As  long 
the  boundary  condition  is  sati  Led  on  the  ship  hull,  the 
calculation  can  still  be  done. 

eik?   =  eikx  [1  +  (eik(^x)  -  1)]  (50) 

From  equation  (43),  neglecting  the  time  dependence  and 
taking  the  factor  e  '^'   outside  the  integral  sign,  it  results: 
ikx    m 

1  n=o   L   n     n 

+  /  A  (U  [elk(?~x)  -1]  r\>      SO,  (51) 

which  is  an  interesting  form  for  our  analysis.   If  only  the 
i       first  integral  existed,  the  result  would  be  the  same  as  the 

one  obtained  applying  the  usual  strip  method  to  the  heave 
j   motions.   The  second  integral  causes  an  additional  deformation 

of  the  wave  and  constitutes  a  correction  factor. 
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A  first  approach  to  this  required  correction  can  be 
calculated,  and  leads  to  the  conclusion  that  the  oncoming 
wave  is  more  deformed  as  it  progresses  from  the  bow  to  the 
stern,  resulting  in  a  decrease  of  the  effective  wave  height. 

The  wave  profile  using  as  reference  the  plane  z     is 
calculated  from  equation  (43)  using  the  following  i     ionship: 

"V-I'M'-o         .  <52> 

The  values  given  by  (52)  should  therefore  be  used  in  the 
calculation  of  the  loading  function  (34). 

lysii      carried  out  for  a  ship  going  in  head 
seas.   The  derivation  of . the  shape  of  the  distorted  wave  at 
ship  side  is  strictly  valid  for  oj->0  ,  but  acceptable  results 
are  expected  when  using  the  same  analysis  for  w^O. 

For  oblique  wav^s  an  ;  ct:    :  in 

considering  instead  of  equation  (30) 


2tt 
e 


w~  =  j~    ^V  +  c°cos  a)  (53) 


where  A   is  an  effective  wave  length 
e 

X      —   X •  s  e  c  a  (54) 

e 

and  X    is  the  wave  length,  V  the  ship's  speed,  c  the  wave 
celerity  and  a  is  the  angle  between  the  ship's  heading  and 
the  wave  direction,  with  a=0  for  head  seas. 
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2.4-PROCEDURE 

Since  the  mathematical  procedures  of  Sections  2 . 2  and  2 . 3 
are  rather  long,  it  was  felt  it  would  be  convenient  to  present 
in  an  abridged  form  the  results  obtained.   The  purpose  of  this 
section  is  to  provide  such  a  summary,  as  well  as  indicate  how 
the  different  parameters  in  the  equations  can  be  calculated. 

The  bending  moment  due  to  springing,  evaluated  amidships, 
is  given  by  equation  (32)' 

0 

M  =  -co?  /  _  /0  x-y(x)  °X.  (x)  'dx-q.  (t)  (32) 

s      i   -L/2  i        Ji   '  ' 

or  equation  (33) 

0 
Ms  -   /_L/2  x-u(x) •X1(x) •dx-q1(t)  (33) 

The  effective  mass  per  unit  length  \i    is  the  sum  of  the 
lip's  mass  pei  unil  3  agth  (taken  from  the  shir'  : 
distribution)  with  the  added  mass  (high  frequency  limit) . 
The  added  mass  can  be  calculated  using  the  Frank  close-fit 
method  [32] . 

The  normal  mode  shapes  X. (x)  can  be  calculated  using  the 
Prohl-Myklest&d  digital  method  mentioned  in  [33]  or  t) 
finite  element  pre  [34]. 

Equation  (31)  gives  a  suitable  expression  for  q. (t) : 
Q  exp(j   w  t) 

q,(t)  =  — — - - —  (3D 

PC  -to2u  .  +j    to      c. 

i      e    i    J      e      l 

with 

Qi(t)    -   QQ   exp(j      to      t)  (29) 


and 

L/2 
Qi(t)  =  f_L/2    P(x,t)'X±(x)  dx  (11 


For  the  frequency  of  encounter  to   expressions  (30)  or 

(53)  can  be  used.  K.  is  defined  by  equation  (24): 

K.  -  to2  y.  (24) 

ill 

where  the  natural  frequencies  to .  can  be  computed  using 

reference  [34]  or  the  Prohl-Myklestad  method  reported  in  [33]; 

the  deiinition  of  y..  is 


L/2 
U±  =   /_L/2   M(x).X?(x)   .  (16) 


C.  is  given  by  equation  (26)- 
c_  L/2    dA 

c  ,  =  (— )  y.  •  v  /  _  .,  -  .  -    x?  (  )  ■ .:  (;  ) 

i     y     l        l/2  i 


where  V  is  the  ship's  speed  and  A__  is  the  added  mass  (high 
frequency  limit) ,  which  can  be  calculated  using  reference 
[32].   The  structural  damping  coefficient  (c  /y )  can  be 
taken  from  figure  2  of  reference  [10],  of  which  the  following 
exp3  •-   '  c  n  ::   a  close  fit ; 


fi   =  (_^_)9/5  ,    . 

y      l  50  ;  (       ' 

The  loading  function  P(x,t)  is  given  by  equation  (34); 

n  2iTh 

dA    (  — j— 

P(x,t)=-CpgBn*+[N  (x)-V~^-]n*+A_3n*}e  (34) 
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Here  p,  g,  B,  V  and  A  are  known.   The  local  damping  force 

coefficient  N  (x)  and  the  local  sectional  added  mass  A__  can 
z  33 

be  evaluated  using  reference  [32].   The  mean  draft  h  is 
approximated  by 


h  -  H  Cs  (35) 

with  H  the  local  section  draft  and  C  the  local  section 

s 

area  coefficient. 

The  effective  wave  profile  n*  is  calculated  from  (52) : 

n*  =  ~&y:l„__n  (52) 

g  dt  z— 0 


.and  (43)  : 

exp(io)  t)    <»  L/2 
<f»(x,y,z)  = 5 — -     2  /   .  Ate)  ^n[(e-x)  ,yfz]d£    (43) 

n.   Q 


The  constants  A  (£}  are  calculated  from  equation  (38)  u: 

the  boundary  condition  (42);  t|)   is  defined  in  equations  (44) 

and  (45). 

In  (34),  the  derivatives  n*  and  n*  are  calculated  using 

defdr  ;  irons  (36)  and  (37): 

Dn  * 


and 


Dn*     3      9 

n*  =  ^-  =  |£-  v^.)n*  (37) 
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3.   CONCLUDING  REMARKS 

3.1-DESIGN  CONSIDERATIONS 

Taking  a  closer  look  at  the  mathematical  model  developed 
above,  it  is  possible  to  draw  some  conclusions  about  the 
importance  of  some  design  parameters  on  the  occurrence  and 
magnitude  of  springing. 

Equation  (33),  which  gives  the  midships  bending  moment 
due  to  springing,  shows  that  a  short  hull  is  better  than  a 
long  one.   This  is  reflected  both  on  the  limits  of  integration 
and  on  the  mode  shape  X, (x) ,  which  tends  to  be  larger 
(amids]  Lps  and  at  both  snds  of  the  hull)  for  a  long  ship. 
The  same  conclusion  can  be  drawn  using  the  same  reasoning 
from  equation  (31)  when  (29)  and  (11)  are  substituted;  note 
that,  for  a  sma"1!  q.  (t)  ,    a  high  C.  is  desirable,  and  equation 
(26)  indicates  again  the  advantage  of  a  short  hull. 

The  same  conclusion  about  bhe  ship  length  is  obtained  if 
the  hull  natural  frequency  of  vibration  is  considered. 

The  empirical  equation  (56)  is  a  well-known  form  for  the 
two-mode  natural  frequency  of  ships'  hulls  [35]: 


where : 

N  -  frequency  per  minute 

I  -  moment  of  inertia,  in  iri2ftz 

A  -  ship  displacement  in  tons 

L  -  ship  length  in  feet 

(j,  -  empirical  coefficient 
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Therefore  the  natural  frequency  is  increased  by  increasing 
I  (by  increasing  the  ship's  depth  or  using  mild  steel),  or 
by  decreasing  the  product  AL 3  at  a  higher  rate  than  I.   Since 
at  higher  frequencies  the  wave  energy  spectrum  has  smaller 
values,  the  excitation  will  be  smaller.   This  explains  the 
higher  values  of  springing  on  ships  with  high  block  coeffici- 
ents and  long  hulls;  an  important  increase  in  springing  is 
expected  in  ships  longer  than  900  ft  LBP  [10] . 

Considering  y(x),' which  is  the  sum  of  the  ship's  m?ss 
with  the  added  mass  (per  unit  length) f  in  equation  (33)  it  is 
conclude       a  r      u(x)  is  desirable  sr  le   Lly   .  the 
points  of  the  hull  wh   -  Xn  (x)  is  larger.   Therefore  a  fairly 
uniform  weight  distribution  leads  to  situations  more  suscep- 
tible to  higher  springing  than  a  ship  with  the  weight  amid- 
ships (if  the  weigh  mass  could  be  concentrated  at 

of  vibration,  then  there  would  be  no  spri] 
bending  moment  amidships . ) 

It  is  interesting  to  note  here  that  calculations  carried 
out  assuming  a  uniform  weight  distribution  along  a  ship 
length  lead  to  a  natural  frequency  24%  lower  than  the  one 
found  using  the  ■  .  istribution  [36] .   Since 

lower  natural  frequencies  give  higher  springing  moments, 
these  calculations  reinforce  the  conclusion  mentioned  on  the 
last  paragraph. 

It  was  said  before,  and  noticed  on  equations  (33)  and  (26) 
that  a  small  added  mass  A_~  is  desirable.  Experiences  carried 
out  with  beams  of  rectangular  cross  sections  sho;v  that  the 
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added  mass  increases  as  the  ratio  beam/depth  increases.   For 
low  values  of  springing,  it  is  therefore  concluded  that  a 
small  beam/depth  ratio  is  an  advantage. 

Summarizing  the  considerations  discussed  above,  the 
following  general  recommendations  can  be  made  in  order  to 
obtain  a  design  that  will  experience  smaller  springing 
stresses : 

— short  hulls 

— small  block  coefficient 

— weight  distribution  concentrated  in  the  cent   1  part 
o  t]  s  hu]  1 

--small  beam/depth  ratio 

The  relative  importance  of  the  different  parameters 
mentioned  cannot  be  assessed  at  this  point,  and  should  be 
investigated  by  computing  the  springing  response  for  various 

ign  cases  and  using  model  tests  as  well  as  full  scale 
experiments . 

3 . 2-METHOD  EVALUATION 

It  is  important  to  note  that  the  springing  phenomenon 
from  waves  which  are  she  t  relative  to  the  ship 
length,  while  the  theorectical  considerations  to  evaluate 
the  wave  generated  forces  on  a  ship  were  derived  for  conditions 
of  long  wave  lengths  compared  to  the  dimensions  of  the  ship 
cross-section.   Therefore  the  question  of  the  validity  of 
such  calculations  remain  until  some  further  research  is  done 
on  this  matter,  because  the  analytical  computations  are 
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dependent  upon  the  accuracy  of  theorectical  expressions  that 
have  not  been  completely  verified  for  the  range  of  high 
frequency  at  which  springing  occurs. 

The  mathematical  model  that  is  the  subject  of  the  present 
study  should  be  evaluated  by  applying  it  to  different  design 
cases  for  which  experimental  data  is  available.   It  is 
suggested  that  a  computer  program  be  written  in  order  to  allow 
the  calculation  of  the  springing  response  of  different  vessels 
to  different  excitations.   As  an  example,  it  is  recommended 
that  t       od  be  applied  to  the  cases  :  ported  in  [37],  where 

y  important  experin         ba  c  •    s\   ringing  is  collecti  3 
and  discussed. 

Future  investigations  on  springing  response  should  use 
'  dir   t  measure  of  the  exciting  waves,  such  as  a  wave 
mounted  over  the    w-  of  the  ship.   These 
would  then  provide  a  set  of  wave  records'  that  would 
directly  correlated  with  the  measured  stresses,  and  would 
lead  to  a  better-  understanding  of  the  exciting  mechanism  of 
the  springing  phenomenon  as  well  as  allow  a  positive 
evaluation  of  the  mathematical  models  proposed  to  calculate 
■  springing  res 
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